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Abstract 

<d ; 

r~| | The top quark has a large Yukawa coupling with the Higgs boson. In the usual 

extensions of the standard model the Higgs sector includes extra scalars, which also 
tend to couple strongly with the top quark. Unlike the Higgs, these fields have a natural 
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, mass above 2mt, so they could introduce anomalies in tt production at the LHC. We 

IT) , _ 

T^j- , study their effect on the tt invariant mass distribution at yfs = 7 TeV. We focus on 

the bosons (H ,A) of the minimal SUSY model and on the scalar field (r) associated to 
the new scale / in Little Higgs (LH) models. We show that in all cases the interference 
with the standard amplitude dominates over the narrow-width contribution. As a 
consequence, the mass difference between H and A or the contribution of an extra 
T-quark loop in LH models become important effects in order to determine if these 
fields are observable there. We find that a 1 fb" 1 luminosity could probe the region 
tan ,3 < 3 of SUSY and v/(y/2f) > 0.3 in LH models. 
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1 Introduction 



The main objective of the LHC is to reveal the nature of the mechanism breaking the 
electroweak (EW) symmetry. This requires not only a determination of the Higgs mass and 
couplings, but also a search for additional particles that may be related to new dynamics 
or symmetries present at the TeV scale. The top-quark sector appears then as a promising 
place to start the search, as it is there where the EW symmetry is broken the most (it 
contains the heaviest fermion). Generically, the large top-quark Yukawa coupling with the 
Higgs boson (h) also implies large couplings with the extra physics. For example, in SUSY 
extensions h comes together with neutral scalar (if) and pseudoscalar (A) fields pQ. Or 
in Little Higgs (LH) models, a global symmetry in the Higgs and the top-quark sectors 
introduces a scalar singlet and an extra T quark [21 13] . In all cases these scalar fields have 
large Yukawa couplings that could imply a sizeable production rate in hadron collisions and 
a dominant decay channel into ti. 

The energy and the luminosity to be achieved at the LHC make this collider a top-quark 
factory, with around 1.5 x 10 5 pairs at -^/s = 7 TeV and 1 fb _1 . In this paper we study 
the possibility that the production and decay of extra Higgses distorts the ti invariant mass 
distribution (rritf). The relevant amplitudes are pictured in Fig. 1. We first review (U E] 
the (analytical) expressions for the cross section when the intermediate field is a scalar or 
a pseudoscalar field and the loop fermion is the top or a heavier T quark. Then we define 
the models to be analyzed and study the parton-level cross section in each case. Finally we 
discuss the possible signal at the LHC. 

2 Top quarks from scalar Higgs bosons 

The potential to observe new physics in m t i at hadron colliders has been discussed in previous 
literature [5J El [7J EJ |9j [101 111 Il2l EE]- I n general, any heavy s-channel resonance with a 
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Figure 1: Diagrams that interfere in tt production. 
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significant branching ratio to ti will introduce distortions: a bump that can be evaluated in 
the narrow-width approximation or more complex structures (a peak followed by a dip) when 
interference effects are important [14]. In the diagram depicted in Fig. 1 the intermediate 
scalar is produced at one loop [15] , but the gauge and Yukawa couplings are all strong. 

Let us first consider a scalar coupled to the top quark and (possibly) to a heavy fermion. 
The leading-order (LO) differential cross section for gg — > ti is then 
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where z = cos 9 is the cosine of the angle between an incoming g and t, m t and y^it are the 
top-quark mass and Yukawa coupling, and (3 - 



'1 — Am 2 / s is the velocity of t in the center 
of mass frame. The function N{s) associated to the fermion loop is 
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where / may be the top or another quark strongly coupled to 
form depending on the mass mj: 
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If 2m f > yfs then If is real and the interference vanishes at s = m 2 ^. If / is the top or any 
fermion with 2m f < y/s, then this contribution can be seen as a final-state // interaction 
[5]. The diferential QCD contribution daQco/dz can be found in [i~6l [T7] . 

For a pseudoscalar A we have 
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Figure 2: a(gg — > tt) with a standard Higgs of mass m h = 500 GeV. 



As we will see in the next section, to have an observable effect it is essential that the 
width is small. This is precisely the reason why the effect on m t i of a very heavy standard 
Higgs h would be irrelevant. A 500 GeV Higgs boson would couple strongly to the top quark, 
but even stronger to itself: A = m\j (2v 2 ) ~ 2. Its decay into would-be Goldstone bosons 
(eaten by the massive W and Z) would then dominate, implying a total decay width 
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and we have taken a common W, Z mass my f» 90 GeV. 

The plot in Fig. 2 shows a too small deviation due to the standard Higgs in cr(gg — > tt). 
To have a smaller width and a larger effect the mass of the resonance must not be EW. In 
particular, SUSY or LH models provide a new scale and massive Higgses with no need for 
large scalar self-couplings. 
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3 SUSY neutral bosons 



SUSY incorporates two Higgs doublets, and after EW symmetry breaking there are two 
neutral bosons (H and A) in addition to the light Higgs. The mass of these two fields is not 
EW (it comes from the SUSY breaking sector), so they are naturally heavy enough to decay 
in ti. Their mass difference is EW, of order M^/m A , with important top-quark corrections 
at the loop level. More precisely, the relation between m A and is [TS] 
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Varying the /i parameter and the stop masses and trilinears, for m A = 500 GeV we obtain 
typical values of m# — m A between —2 and +10 GeV. 

The scalar masses of interest correspond to the decoupling regime, where h is basically 
the SM Higgs and 

m t V2 1 

Vmt ~ 7 a w • (I 6 ) 

v tan p 

In addition, we will consider low values of tan /3, where the decay into bottom quarks is not 
important and the (energy-dependent) widths can be approximated to 
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16tt mn 



IQir m A 



(17) 



22 
21 
20 
19 
b 18 
17 
16 
15 



P- 



■ 1 

1 
1 







i : 
\ : 


- / / . 




/ 




// 




//' 




// 




r 




: // 




■ 


/ 





400 450 500 550 
T (GeV) 



600 




400 450 500 550 
(GeV) 



600 



Figure 3: <j(gg — > tt) for tan/3 = 2 and SUSY bosons of mass itia = fnH — 500 GeV (left) 
or rriA = 500, mn = 505 GeV (right). Dashes provide the narrow- width approximation and 
dots the standard model cross section. 




Figure 4: Standard angular distribution for the t quarks from qq and gg collisions at = 500 
GeV. We include (dashes) the distribution from gg at the peak and the dip of Fig. 3 left . 
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In Fig. 3 we plot cr(gg — > tt) at center of mass energies around rriA = 500 GeV for 
mil = 500 GeV (left) and m# = 505 GeV (right). We have taken tan/3 = 2, which implies 
T H 3.0 GeV and Ta ~ 5.3 GeV. We observe an average 5.5% excess and 8.1% deficit in the 
5 GeV intervals before and after y/s = 500 GeV, respectively. We include in dashes the result 
ignoring the interference (last term in Eqs. fl2]) and (|6])), which would not be captured if one 
uses the narrow-width approximation. It is apparent that the interference with the standard 
amplitude gives the dominant effect. In Fig. 3-left the position of the peaks and dips caused 
by H and A overlap constructively (notice, however, that in this CP conserving Higgs sector 
their amplitudes do not interfere). In contrast, in Fig. 3-right their mass difference implies 
a partial cancellation between the dip caused by A and the peak of H. 

The scalar and pseudoscalar couplings with the top quark grow at smaller values of tan j3, 
increasing the cross section and the scalar width. For example, for tan = 1 the excess at 
y/s < 500 GeV grows to the 6.2% and the deficit to the 9.7%, whereas for tan/3 = 5 the 
excess and deficit are just a 2.1% and a 2.6%, respectively. 

The normalized angular distribution of the t quark in the center of mass frame is given in 
Fig. 4. We plot the standard distributions for top-quark production in gg and qq collisions 
together with the distribution from gg at the peak and the dip obtained in Fig. 3-left. In 
the narrow-width approximation a scalar resonance gives a flat contribution. However, we 
find that the excess or deficit from the scalar interference is not flat and does not change 

I I I I 18.0 i 
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Figure 5: c(gg — > ti) for uia = 400 GeV and rriA = 700 GeV. We include the cases 
m H = 400,408 GeV (left) and m H = 700, 703 GeV (right) 
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significantly the angular distribution. Different cuts could be applied to reduce the back- 
ground for tt production at the LHC [7] or even to optimize the contribution from gg versus 
qq, but not to enhance the relative effect of the scalars on a(gg —> tt). 

In Fig. 5 we plot the parton cross section for lower and higher values of the pseudoscalar 
mass (rriA = 400,700 GeV). We include the cases where the boson H is degenerate with 
A or slightly heavier (m# = 408 GeV and m# = 703 GeV). We see that at lower scalar 
masses the peak dominates, whereas for large values of tua the dip is the dominant effect. 
This behaviour, related to the slope of the standard cross section, reduces in both cases the 
relevance of the mass difference between the scalar and the pseudoscalar Higgses. 



4 Little Higgs boson 

In LH models the Higgs appears as a pseudo-Goldstone boson of a global symmetry broken 
spontaneously at the scale / > v/y2 = 17 A GeV. The global symmetry introduces an 
extra T quark that cancels top-quark quadratic corrections to the Higgs mass parameter. 
The presence of this vectorlike T quark and of a massive scalar singlet (the Higgs of the 
symmetry broken at /) are then generic features in all these models. 

Once the electroweak VEV is included the doublet and singlet Higgses (and also the t 
and T quarks) mix [191 [20]. The singlet component ~ v/(y2f) in h will reduce its coupling 
both to the top quark and to the gauge bosons and, in turn, r will get a doublet component 
that couples to these fields. 

It is easy to see that the most general] top-quark Yukawa sector with no quadratic 
corrections at one loop is 
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lr Thcrc could be an additional mixing, T — > cpT + spt in the second line of Eq. (|18[) . but it must be small 
[21] to avoid a too large value of Vrb ■ 
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Figure 6: cr(gg — >■ £i) for a LH model with se = 0.5, m r = 500 GeV and uit = 400, 500, 700 
GeV. Dashes (dot-dashes) correspond to an amplitude with only the T (t) quark loop. 

A (/ + ^) (c cos - s e si, TT< + h.c. . (20) 

Fermion masses, Yukawa couplings and dimension-5 operators (necessary to check the can- 
cellation of all one-loop quadratic corrections) are then obtained by expanding 

h h 2 h h . . 

cos—^Ril-—T, sin-^w-p- (21) 
y/2f 4/ 2 x/2/ V2f 1 ' 



The fermion masses and the Yukawas to the heavier scalar r have the same structure, 
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/ v f 

where the quarks are mass eigenstates. The mass of the heavier T quark is m T m t cg/(sgc a ), 

and its mixing with the doublet Vrb ~ s 2 e s a c a /c^. 

The extra Higgs r is somehow similar to the heavier scalar in a doublet plus singlet model, 
with the doublet component growing with sq = v/(\/2f). If sg is sizeable so is its coupling 
to the top quark. The coupling to the extra T quark is stronger, but if r is lighter than 
2rriT then its main decay mode will be into tt. Actually, the doublet component in r may 
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also imply large couplings to the would-be Goldstone bosons for large values of m r . More 
precisely, its decay width T r (s) at 4m^ < s < 4m^ is 



Therefore, r is a naturally heavy (m r m f) but narrow scalar resonance with large couplings 
to quarks and an order one braching ratio to ti. 

In Fig. 6 we plot the parton cross section a(gg — > ti) for se = 0.5, ttit = 500 GeV and 
several values of m r . We separate the contributions from the top and the T quark loops (the 
second one vanishes at s = m%). The plot is similar to the one obtained for SUSY bosons 
of the same mass. At higher values of m r the decay width T r grows, diluting the effect (see 
Fig. 6, right). In contrast, for lower masses the scalar r has a narrow width and is strongly 
coupled to quarks, which produces a larger effect (in Fig. 6, left). The contribution from the 
standard t-quark loop grows with sg, whereas the contribution from the extra T-quark is 
basically independent of mr- 



Let us now estimate the invariant mass distribution of ti events (mu) in pp collisions at the 
LHC. To evaluate the hadronic cross sections we will use the MSTW2008 PDFs [22]. The 
effect of next-to-leading order (NLO) corrections on the expressions given in previous sections 
has been studied by several groups (see for example [6j [23] ) . In particular, the authors in [6] 
analyze the dependence of da/dm t t on the choice of renormalization and factorization scales 
and of PDFs. They show that if the LO cross section is normalized to the NLO one at low 
values of m t t, then the deviations introduced by these scales and by the uncertainty in the 
PDFs at m t i < 1 TeV are small (order 10%). For (scalar and pseudoscalar) Higgs production 
in pp collisions and Higgs decay, a complete review of NLO results can be found in [1] . From 
the expressions there we obtain that QCD corrections enhance the production cross section 
in approximately a 20%, and that the Higgs decay width into ti (for 3> 2m t ) is also 
increased in around a 10%. Given these estimates, we have evaluated pp — > ti taking fixed 
renormalization and factorization scales (/j-r,f — m t) an d normalizing the LO result to the 
NLO cross section in [6] with a global factor of 1.3. Our differential cross section coincides 
then with that NLO result at m t j = 500 GeV. 

We will take a center of mass energy of 7 TeV. We obtain that at these energies the 
cross section pp — > ti is dominated by gg fusion, with qq — > ti accounting for just 12% of 
the top-quark pairs. In Figs. 7,8 we plot da/dm t t for some of the SUSY and LH models 




+ 



^(^ + ^(1-^) 2 ) 



(24) 



5 Signal at the LHC 
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Figure 7: da/dm t i in SUSY models with m H = rriA = 400,500,700 GeV and tan/3 = 2 
(solid), 5 (dashes). 
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Figure 8: da/dm t t in LH models with m r = 400,500,700 GeV and s e = 0.5 (solid), 0.2 
(dashes). 
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Figure 9: Number of tt events in pp collisions at 7 TeV and 1 fb 1 for itla = rriH = 500 GeV 
and tan (3 = 2 distributed in 5 GeV bins. 
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Figure 10: Deviation A = (N—Nsm) / V Nsm i n the number of events respect to the standard 
prediction for two different binning (m^ = itlh = 500 GeV and tan/3 = 2). 
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described in Sections 3-4. These figures translate the parton-level cross sections in Figs. 3-6 
into anomalies in the invariant mass distribution in pp collisions. 

To estimate the possible relevance at the LHC of these cross sections, we will calculate 
the number of ti events assuming an integrated luminosity of 1 fb _1 (expected at the 7 TeV 
phase of the collider), and we will not apply any cuts. In Fig. 9 we plot the number of events 
per 5 GeV bin of m t i in the SUSY model with = m H = 500 GeV and tan/3 = 2. We 
observe a 5% excess followed by a 9% deficit, with smaller deviations as m t i separates from 
the mass of the extra Higgs bosons. In Fig. 10 we distribute the events in 20 GeV bins and 
plot the statistical significance 

of the deviations, where N is the total number of events in the bin. The typical signal is 
an increasing excess in a couple of 20 GeV bins that may reach a +3. 4a deviation followed 
by a deficit of —4.5a. We find that changing the binning is important in order to optimize 
the effect. If the same 20 GeV bin includes the peak and the dip (Fig. 10, right) then the 
maximum deviation is just a ±2a effect. 

The result is very similar for a LH scalar of m r = 500 GeV with sg = 0.5. In this LH 
model we obtain deviations in consecutive 20 GeV bins reaching +2. 5a and — 2a. However, 
the effect is a bit more localized, and the cancellation if peak and dip coincide in a bin is 
stronger: it may result in three bins with just +1.3(7, +0.6a and — 1.2a deviations. 

The binning is less important for larger Higgs masses. For example, in the SUSY case 
with rriA = inn = 700 GeV the typical sequence is a couple of 20 GeV bins with a slight 
+0.2a excess followed by —1.2a, —0.4a and —0.2a deficits. In the LH model with m r = 700 
GeV the initial excess (caused by the T-quark loop) is a bit more significant, a typical 
sequence would consist of two bins with +0.4a excess followed by —0.8a and —0.4a deficits. 

Let us finally focus on lighter Higgses, as they provide the most promising signal. In 
Fig. 11 we plot the event distribution (left) and the statistical significance (right) for tan/3 = 
2 and = itlh = 400 GeV, whereas Fig. 12 corresponds to a mass difference of 8 GeV 
(itla = 400 GeV and van — 408 GeV). The sequence of deviations in both cases would be 
seen as a clear anomaly, reaching an excess of up to 13a (for itlh — ttla = —2 GeV) in a 
20 GeV bin. The LH case is analogous but, again, more localized. We obtain an excess of 
+3. 4a in a 20 GeV bin followed by a —1.7a deficit. 
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Figure 11: Number of tt events in pp collisions (left) and deviation A (right) for wla = itlh 
400 GeV and tan (3 = 2. 
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Figure 12: Number of tt events in pp collisions (left) and deviation A (right) for itla = 400 
GeV, m H = 408 GeV and tan/3 = 2. 
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6 Summary and discussion 



In models with an extended Higgs sector the extra bosons tend to have large couplings 
with the top quark that imply a sizeable one-loop production rate at hadron colliders. If 
the mass of these bosons is not EW but comes from a new scale (e.g., the SUSY or the 
global symmetry-breaking scales), then they may decay predominantly into tt. We have 
studied their effect on the tt invariant mass distribution at 7 TeV and 1 fb _1 . We have 
considered the deviations due to the neutral bosons A and H of the MSSM, and to the 
scalar r associated to the scale / in LH models. In all cases the interference dominates, 
invalidating the narrow-width approximation. 

The effect for masses around 500 GeV is a peak followed by a dip of similar size. In 
the SUSY case, values of mn — ttia smaller than 3 GeV enhance the deviation, whereas 
for larger values the effects at < m t t < inn tend to cancell each other. For tan (3 = 2 
the significance of the signal, that can be optimized by changing the binning, results in 
sequences of 20 GeV bins with +1.4(7, +3.4<r, — 4.5<r, —1.0a deviations. In LH models the 
field r couples both to the top and to an extra T quark. The main difference with the SUSY 
case is that the T quark is heavy and the one-loop form factor to produce the scalar does not 
get an imaginary part. To have a significant effect the doublet component sg — vj (V2f) in r 
must be large and the extra T quark heavier than m r /2 (to close the TT decay channel). As 
sg grows r resembles the standard model Higgs, but with a singlet component that reduces 
its width. The signal at the LHC for sg = 0.5 and m r = 500 GeV is similar to the SUSY 
case just described. 

At larger scalar masses the peak decreases and the effect is basically a dip in the invariant 
mass distribution. For = 700 GeV and tan (3 = 2 we get a couple of 20 GeV bins with a 
deficit of — 1.2<7, and — 0.4o\ The effect that one may expect in LH models is alike, although 
(due to the contribution of the T-quark loop) the difference between peak and dip is smaller. 

Lower scalar masses provide the most promising scenario. Here the peak dominates both 
in SUSY and LH models. In the SUSY case with itla = 400 GeV and uih = 408 GeV the 
sequence of 20 GeV bins at 7 TeV and 1 pb _1 consists of +2. la, +9.9(7, — 4.2cr and — 2.3<r 
deviations. The signal increases in up to a 30% if the mass difference is smaller (the optimal 
case is uih —iTiA = —2 GeV). In the LH case with sg = 0.5 and m r = 400 GeV the sequence 
reaches a +3. 4a deviation. All the effects grow for lower values of tan (3 and of the LH 
scale /. Obviously, their observability will be better if the LHC reaches 14 TeV and higher 
luminosities. 

An important observation is that the angular distribution of the t quark is unaffected 
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by the intermediate scalar or pseudoscalar bosons. The excess or the deficit caused by its 
interference with the standard amplitude does not have a flat distribution in the center of 
mass frame, as one obtains in the narrow-width approximation. 

Finally, we would like to comment on the possibility to observe this type of signal at the 
Tevatron, which may achieve 10 fb _1 at ^/s = 1.96 TeV. The main difference with the LHC is 
that at the Tevatron 90% of the top-quark pairs are produced through qq interactions. Since 
the signal that we have explored is caused by interference in the gg — > tt channel, for the 
same integrated luminosity the deviations there would be 9 times weaker than at the LHC 
(where gluon fusion provides 90% of the top pairs). We find, however, that la deviations 
could be obtained at the Tevatron for low masses of the heavy Higgs bosons. This signal 
could be enhanced by separating the tt events in two or three sets according to the cos 9 of the 
final t quark. As we see in Fig. 4, the gg and qq contributions at m t i 500 GeV have 

different angular distributions (this difference, however, vanishes at lower invariant masses). 
One could separate, for example, the events with |cos 6 1 ! larger or smaller than 0.6. Then the 
anomalies in da /dm t i that we have discussed should increase in the |cos#| > 0.6 interval. 

Although the generic effect on the tt invariant mass distribution caused by a scalar field 
with strong couplings to the top quark is known, we think that it is also important to study 
particular models. We find that the peculiar anomalies that appear in these two extensions 
of the standard model may be optimized by changing the binning and by applying the same 
cuts that select ti production from gluon fusion. 
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